DFTT 7/11 
HIP-2011-11/TH 



Thermodynamics of SU(N) Yang-Mills theories in 2 + 1 
dimensions I — The confining phase 



Michele Caselle a , Luca Castagnim b , Alessandra Feo a , Ferdinando Gliozzi a and Marco Panero c 



Dipartimento di Fisica Teorica dell'Universita di Torino and INFN, Sezione di Torino, 

Via P. Giuria 1, 1-10125 Torino, Italy 

b Institute for Theoretical Physics, University of Regensburg, 
D-93040 Regensburg, Germany 

c Department of Physics and Helsinki Institute of Physics, University of Helsinki, 

FIN-00014 Helsinki, Finland 

E-mail: caselleOto . inf n. it , luca. castagnini@physik.uni-regensburg.de , 
feo@to.infn.it, gliozzi@to.infn.it, marco.panero@helsinki.fi 



Abstract 

We compute the equation of state in the confining phase of SU(iV) Yang-Mills theories 
with N = 2, 3, 4, 5 and 6 colors in 2 + 1 dimensions, via lattice simulations. At low enough 
temperatures, the results are accurately described by a gas of glucballs, including all known 
states below the two-particle threshold. Close to the deconfinement temperature, however, this 
prediction underestimates the numerical results, and the contribution from heavier glueballs 
has to be included. We show that the spectral density of the latter can be accurately described 
using a simple bosonic string model. 



PACS numbers: ll.10.Wx, 11.15.Ha, 12. 38. Aw, 12.38.Gc, 12.38.Mh 



1 Introduction and motivation 



Determining the phase diagram of strongly interacting matter is a major challenge in elementary 
particle physics, both theoretically and experimentally. 

From the theoretical point of view, the qualitative expectation that, when the temperature or 
the density is sufficiently high, usual hadronic matter gives way to a state of deconfined particles, 
is a straightforward consequence of asymptotic freedom, and has been around since the early days 
of QCD pQ. However, a derivation of the quantitative details of the QCD phase diagram is com- 
plicated by the fact that perturbative methods in thermal gauge theories are typically hindered by 
severe infrared divergences [2], and cannot be reliably applied close to the deconfinement point, 
where the physical coupling is not small. This leaves numerical computations based on the lattice 
regularization of QCD as the main tool for a first-principle study of the QCD phase structure as 
a function of the temperature T, and for vanishing or small values of the quark chemical potential 
[i: recent results are summarized in refs. [3]. 

On the other hand, on the experimental side, the creation of a deconfined plasma of quarks 
and gluons (QGP) in the laboratory has been the goal of a three-decade-long programme of 
heavy ion collision experiments, first at AGS and SPS, then at RHIC, and currently at LHC. In 
particular, the SPS, RHIC and LHC runs have provided convincing evidence for the creation of 
a state of deconfined matter, which achieves rapid thermalization, and can be characterized as a 
nearly ideal fluid [3J. 

It is important to point out that, in these experiments, the physical features of the decon- 
fined plasma are studied indirectly, namely, they are reconstructed from the properties (yields, 
momentum distributions et c.) of the hadrons produced after the expansion and freeze-out of 
the "fireball". The statistical analysis of these results is based on the assumption that, below 
the characteristic temperature range where deconfinement takes place (approximately between 
150 and 190 MeV), thermal QCD can be modelled as a gas of massive, non- interacting hadronic 
resonances [5]- In fact, the very idea of a deconfinement temperature, as the limiting upper 
temperature at which the exponential growth of the density of states in the hadronic spectrum 
would lead to a divergence of the partition function, is even older than QCD [6] . 

Since hadrons are intrinsically non-perturbative objects, any first-principle test of the thermo- 
dynamic description of the confining QCD phase via the hadron resonance gas model necessarily 
requires lattice simulations. It should be noted that this is a computationally challenging task, 
because in the confining phase all equilibrium thermodynamic observables (such as pressure, en- 
ergy density and entropy density) take values, which are much smaller than in the deconfined 
phase. However, the steady increase in computer power and major algorithmic improvements 
have now driven lattice QCD into an era of precision calculations, making it possible to reliably 
investigate the fine details of physical observables, even with limited computational resources. 

Moreover, the most demanding technical aspects in lattice QCD computations involve the 
regularization of fermions, and thus can be easily bypassed, by restricting one's attention to 
the pure-glue sector — which captures most of the features of the full theory, at least at the 
qualitative or semi-quantitative level. The hadronic spectrum of pure Yang-Mills theories has 
been investigated extensively in highly accurate lattice computations [3(8], and several low- lying 
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states are by now well-known: in particular, the lightest state in the SU(3) spectrum is a glueball 
with quantum numbers J PC = ++ and mass (in physical units) about 1.4 GeV — significantly 
heavier than the lightest mesons in the physical QCD spectrum. The masses of glueballs with 
different quantum numbers are also known, and the most recent lattice calculations provide 
precise results for several excited states, too\j 

Looking at pure Yang-Mills theories also offers the further advantage of a cleaner symmetry 
pattern: as opposed to the full-QCD setup, global transformations associated with the center of 
the gauge group are an exact symmetry of the Lagrangian, whose spontaneous breakdown can 
be studied by looking at the expectation value of the associated order parameter, namely the 
trace of the Polyakov loop. The latter provides an unambiguous definition of the deconfinement 
temperature T c , which separates the center-symmetric, confining phase at low temperatures, from 
the deconfined phase at high temperatures. 

The restriction to the Yang-Mills sector is also relevant in the 't Hooft limit of QCD, i.e. in 
the double limit when the number of colors N tends to infinity, and the coupling g 2 tends to zero, 
with the 't Hooft coupling A = g 2 N fixed |10j . In this limit, elementary combinatorics arguments 
show that quark-dynamics effects are subleading (more precisely: suppressed by powers of iV -1 ) 
with respect to contributions involving gluons only, and generic amplitudes for physical processes 
can be rearranged in double series, in powers of the 't Hooft coupling A, and of iV — revealing 
a striking similarity to an analogous expansion in closed string theory: see, e.g., ref. [11] for a 
discussion. While these observations date back to more than thirty years ago, it is interesting 
to note that the large- N limit also plays a technically important role in more modern analytical 
approaches to strongly coupled systems, based on the conjectured correspondence between gauge 
and string theories [12]: according to this correspondence, the string-theoretical dual of a gauge 
theory simplifies to its classical gravity limit, when both the 't Hooft coupling and the number 
of colors in the gauge theory are taken to be large. 

Finally, the emergence of a Hagedorn-like spectrum (i.e. an exponential growth in the number 
of hadronic states, as a function of their mass) has been studied in the large- N limit of QCD, in 
both D = 2 + 1 and 3 + 1 spacetime dimensions, in a very recent work |13j . 

With these motivations, in this paper we report our investigation of the equation of state in 
the confining phase of Yang-Mills theories in 2 + 1 dimensions: we compare our results with a 
hadron resonance gas, using the glueball masses directly extracted from lattice simulations [14j . 
as well as a bosonic string model for the hadronic spectrum. This also allows one to achieve a 
better understanding of the many non-trivial features of effective string models — see, e.g., ref. |15j 
and references therein. 

Our computations can be compared with those reported in ref. [16] for SU(3) in D = 3 + 1: 
in fact, our work can be seen as an extension (in the lower-dimensional case) of the latter study, 
to theories with a different number of colors. This is partially motivated by recent works [17j . 
revealing that the equation of state of the deconfined gluon plasma is characterized by a very 
mild dependence on the number of colors — up to a (trivial) proportionality to the number of 

1 While the extraction of the latter involves a rapidly increasing computational complexity, the spectral density 
of glueball states at higher energies is expected to be approximately described by effective models [9], based on the 
picture of glueballs as closed rings of glue. 
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gluon degrees of freedom — , showing that equilibrium thermodynamic observables in the SU(3) 
theory [18l[T9] are close to the large-iV limit, and lending support to computations based on 
holographic methods and/or on quasiparticle approaches [22]. By contrast, for temper- 

atures T < T c , confinement into color-singlet hadrons leads to the expectation that the number 
of physical states (and the equilibrium thermodynamic quantities) should scale as O(N ), i.e. be 
independent of N, in the large- A r limit. On the other hand, our choice to look at the D = 2 + 1 
setup (rather than D = 3 + 1) is motivated by an important technical aspect: while the decon- 
finement phase transition of SU(3) Yang-Mills theory in 3+ 1 dimensions is of first order, in 2 + 1 
dimensions it is a second-order one. As a consequence, it is expected that the Hagedorn temper- 
ature Tjj should be the same as the deconfinement temperature T c , thereby removing the Th/T c 
parameter to be fitted from the data, and providing a more stringent test of the description of 
the glueball spectral density through a string model. Another lattice study of the equation of 
state in 2 + 1 dimensions (but for the SU(3) gauge theory only) is reported in ref. |23] , 

The structure of this paper is the following: in section [2j we briefly recall the continuum 
formulation and most interesting physical features of SU(A r ) Yang-Mills theories in 2+1 spacetime 
dimensions, define their lattice regularization, and summarize some basic technical information 
about our determination of the thermodynamic quantities on the lattice. In section [31 we present 
the numerical results of our simulations, and compare them with the equation of state predicted 
for a gas of non-interacting glueballs, using the glueball masses known from lattice computations. 
In section HI we define the effective description of the glueball spectrum of SU(iV) Yang-Mills 
theories in 2 + 1 dimensions in the large- N limit through a bosonic string model, and derive 
the corresponding predictions for the equilibrium thermodynamic quantities considered in this 
work. Finally, in section [5] we discuss our findings and their implications. The computation 
of the partition function for an ideal relativistic Bose gas is reviewed in the appendix [XJ while 
appendix [B] reports the derivation of the spectral density for a bosonic closed string model. 
Preliminary results of this study have been presented in ref. |24j . 

2 Non-Abelian gauge theories in 2 + 1 dimensions in the contin- 
uum and on the lattice 

In this section, we first introduce the continuum formulation of SU(iV) Yang-Mills theories in 
2 + 1 dimensions in subsection 12.11 then we discuss their lattice regularization in subsection 12.21 
which also includes some technical details about our computation of thermodynamic quantities. 

2.1 Formulation in the continuum 

Contrary to the D = 1 + 1 case, SU(A r ) gauge theories in D = 2 + 1 spacetime dimensions exhibit 
non-trivial dynamics, and share many qualitative features with Yang-Mills theories in D = 3 + 1. 
They are formally defined through the following Euclidean functional integral: 




(1) 
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In D = 2 + 1 dimensions, the bare square gauge coupling g$ has energy dimension 1, so that 
bare perturbation theory calculations at a momentum scale k are organized as series in powers of 
the dimensionless ratio g$/k |25j . Like in D = 3 + 1 dimensions, also in D = 2 + 1 non-Abelian 
gauge theories are asymptotically free at high energy, and confining, with a finite mass gap and a 
discrete spectrum, at low energy. Their phase diagram as a function of the temperature consists 
of a confined phase (with color-singlet physical states, which can be classified according to the 
irreducible representations of the 0(2) group and charge conjugation) at low temperatures, and 
a deconfined phase at high temperatures. 

The deconfinement transition occurs at a finite critical temperature T c , where the global Z^v 
center symmetry gets spontaneously broken, and the order parameter in the thermodynamic limit 
is the trace of the vacuum average Polyakov loop. In D = 2 + 1 dimensions, the deconfinement 
transition turns out to be a second-order one for SU(2) and SU(3), while it is a very weakly 
first-order one for SU(4), and a stronger first-order one for SU(iV > 5) |26|l27j. 

Equilibrium thermodynamic quantities for SU(iV) Yang-Mills theories in D = 2 + 1 dimen- 
sions can be easily obtained from elementary thermodynamic identities. Let Z(T, V) denote the 
partition function for an isotropic system of two-dimensional "volume" V at temperature T; the 
free energy density /: 

f = -L In Z(T,V), (2) 
is related, in the thermodynamic limit, to the pressure p via: 

p = - lim / . (3) 

V— >co 

In turn, the trace of the energy-momentum tensor A = is related to the pressure by: 

A T'('), (4) 



T 3 dT \T 3 , 

so that the energy and entropy densities (denoted as e and s, respectively) can be expressed as: 

e = A + 2p (5) 

and 

A + 3p 

s = T ■ (6) 



2.2 Lattice regularization 

In this work, we studied non-perturbatively theories based on SU(iV) gauge groups with iV = 2, 3, 
4, 5 and 6 colors, by regularizing them on a finite, isotropic cubic lattice A. In the following, let a 
denote the lattice spacing and l? s x Lt = (N% x Nt)a 3 the lattice volume. The lattice formulation 
regularizes the functional integral in eq. (JT]), trading it for the finite-dimensional multiple integral: 
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where dU a (x) is the Haar measure for each U a (x) G SU(iV) link matrix, and denotes the 
standard Wilson lattice gauge action: 



xgAl<a</3<3 



^ReTrU aP (x) 



2N 

with: = (8) 



where: 

= U a (x)Uf}(x + ofyrfjx + a0)E/J(x). (9) 
Expectation values of gauge-invariant physical observables O are defined by: 

3 

(°) = Y I Uf[au a (x)Oe- s ? (10) 



and can be estimated numerically, via Monte Carlo sampling over a finite set of {U a (x)} config- 
urations; in the following, the number of configurations used in our computations is denoted by 

^conf • 

The numerical results presented in this work are based on sets of configurations (see table [T] for 
details) produced via a Markovian process with local updates; our code implements a combination 
of local heat-bath [28] and overrelaxation steps [29] on SU(2) subgroups [30]. For part of our 
simulations, we also used the Chroma suite [31]. To convert the lattice results obtained from 
simulations into physical quantities, one has to set the physical scale, i.e. to determine the value 
of the spacing a as a function of the bare gauge coupling. This determination is done non- 
perturbatively, via the lattice computation of a reference quantity relevant for low-energy scales 
(such as, for example, the asymptotic slope a of the confining potential V(r) between a pair of 
static sources at zero temperature at large distances r, or the critical deconfinement temperature 
T c ). In this work, the determination of the scale is done using lattice results available in the 
literature [27], and is expressed by the following formula: 

T _ $ - 0.227V 2 + 0.5 

Tc~ Nf (0.357iV 2 + 0.13 - 0.211/iV 2 ) ' ^ ' 

Essentially, the accuracy limits on this formula are set by the precision in the determination of 
the T c ) ' yfo ratio in the continuum limit and in the large- N limit from ref. [27]. The latter reports 
0.9026(23) for the N — > oo limit of this ratio, with a finite-iV correction term proportional to 
iV -2 with coefficient 0.880(43) (this fit is shown to describe well the data, all the way down to 
N = 2). As a consequence, the uncertainty on our determination of the temperature scale can 
be estimated to be of the order of 1%, and has a negligible impact on our analysis (for the sake 
of clarity, we omit the corresponding horizontal errorbars from our plots). 

More generally, it is worth noting that, given that all numerical simulations are done at 
finite values of the spacing a, on the lattice different observables can be affected by different 
discretization artifacts, thus the choice of a particular observable to set the scale introduces a 
systematic uncertainty. However, the quantitative effect of such uncertainty is small, C(a 2 ); for 
a comparison with alternative definitions of the scale, see, e.g., refs. [231 132]. 
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[24.6,32.0] 


1 x 10 5 


6 x 10 5 
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48 2 x 6 


161 


[30.0,46.0] 


2 x 10 4 


1.6 x 10 5 




56 2 x 7 


188 


[34.5,53.2] 


2 x 10 4 


1.6 x 10 5 




64 2 x 8 


200 


[39.0,58.9] 


2.5 x 10 4 


2 x 10 5 
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48 2 x 6 


24 


[51.0,53.2] 


1 x 10 5 


1 x 10 5 








[54.0,58.5] 


1 x 10 5 


5 x 10 5 








[59.0,60.0] 


1 x 10 5 


3 x 10 5 








[61.0,64.0] 


1 x 10 5 


1 x 10 5 


6 


48 2 x 6 


16 


[75.0, 78.0] 


1 x 10 5 


1 x 10 5 








[79.0,83.0] 


1 x 10 5 


4 x 10 5 








[84.0, 86.0] 


1 x 10 5 


2 x 10 5 








[88.0,95.0] 


1 x 10 5 


1 x 10 5 



Table 1: Parameters of the main set of lattice simulations used in this work: N denotes number 
of colors, Nt and N s are, respectively, the lattice sizes along the time-like and space-like directions 
(in units of the lattice spacing), np denotes the number of /3- values (i.e. of temperatures) that 
were simulated, in each f3 m in < /5 < Pmax interval; the T = and finite-T statistics at each 
/3-value are shown in the last two columns. For N > 2, all T = simulations were performed on 
lattices of size (aN s ) s . 
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The lattice simulation of Yang- Mills theories in thermodynamic equilibrium is straightforward: 
the temperature (in natural units h = c = ks = 1) is defined by the inverse of the lattice size, 
T = l/(aN t ), along a compactified direction, with (anti-)periodic boundary conditions for bosonic 
(fermionic) fields, while the sizes in the other directions are kept sufficiently large, N s 3> Nt, to 
enforce a good approximation of the thermodynamic limit — see, e.g., refs. [33] for a discussion. 
To obtain the temperature dependence of all equilibrium thermodynamic quantities (or, more 
precisely, of their difference with respect to the value at T = — for which we run simulations on 
lattices of size (aN s ) s ) we varied the temperature by changing a (which is a function of j3) at fixed 
Nt and N s = 8N t (except for the SU(2) gauge group: see tabled] for details), and repeated the 
computations at increasing values of Nt to estimate discretization effects and perform a continuum 
extrapolation. As compared to the so-called "fixed-scale approach" [34], this method allows one 
to perform efficiently a fine and accurate temperature scan. The trace of the energy-momentum 
tensor can be obtained from: 

A= i^L« u ^-<< u ° M - <i2) 

where (Un)r denotes the expectation value of the average plaquette at the temperature T, while 
the pressure can be determined using the "integral method" [35] : 



p=\f d^((U n ) T -(U n )o), (13) 

a h 



where /3q is a value of the Wilson parameter corresponding to a temperature sufficiently deep in 
the confined phase. In the present work, the numerical evaluation of the integral in eq. (I13p was 
done by the trapezoid rule, except for the SU(4) gauge theory, for which we used the method 
described by eq. (A. 4) in ref. [36], which is characterized by systematic errors 0(n7 4 ). Thus, 
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the uncertainty on the pressure depends on the statistical precision of the plaquette differences, 
and on the systematic uncertainty related to the choice of the lower integration extremum {3q. 
Since the plaquette differences at different values of /3 are obtained from independent simulations, 
the statistical errors on p are obtained using standard error propagation. As for the systematic 
uncertainty related to the choice of the lower integration extremum, we checked that, by virtue 
of the exponential suppression of the plaquette differences in the confined phase, pushing (3q to 
even lower values than those we used, would not have any significant impact on our results for 
the pressure. The energy and entropy densities are then obtained from eq. ([5]) and from eq. Q, 
respectively. 

As a technical remark, note that eq. (j 1 2j) and eq. (|13|) show that the determination of ther- 
modynamic quantities from very fine lattices can be computationally rather demanding, given 
that they are extracted from differences of plaquette expectation values at zero and finite tem- 
perature, and such differences scale like a 3 (or a 4 in the D = 3 + 1 case). However, our numerical 
results reveal a mild cutoff dependence (at least in the range of Nt values that we simulated), 
allowing one to get a reliable extrapolation to the continuum limit. In particular, leading-order 
discretization terms affecting the Wilson lattice gauge action eq. ([8]) are quadratic in a, so that 
continuum results can be obtained by extrapolation of fits in 1/Nf to the Nt — > oo limit. These 
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results can then be compared with the theoretical predictions of effective models for the equation 
of state of Yang-Mills theories in the confined phase, as discussed below. 

3 Numerical results 

Fig. Q] shows our numerical results for the dimensionless ratio of the trace of the energy- momentum 
tensor over the cube of the temperature, A/T 3 , as a function of T/T c , for the various SU(iV) 
gauge groups. In the confining phase, the thermodynamics is that of an ensemble of color-singlet 
states, and, since the number of the latter does not depend on N (with the exception of N = 2, 
for which, due to the (pseudo-)real nature of all irreducible representations of SU(2), there exist 
no states with charge conjugation quantum number C = — 1), it is reasonable to expect that the 
equilibrium observables should not depend strongly on the rank of the gauge group. This is indeed 
observed in fig. [H showing the approximate collapse of data from different groups onto a universal 
curve, up to temperatures around 0.95T C , or even larger. In the same figure, we also show the 
comparison with the curve (dashed line) describing the trace of the energy-momentum tensor 
for a relativistic gas of massive, non-interacting glueballs, using the glueball masses extracted 
from lattice computations in ref. [13], and restricting to states below the two-particle threshold; 
the leading contribution is given by the lightest glueball (dotted line). As it will be discussed 
in section^ both these curves severely underestimate the lattice results at temperatures larger 
than approximately 0.9T C . 

The plot also shows that, very close to the deconfinement transition, the data corresponding 
to different gauge groups start arranging themselves according to the multiplicity given by the 
number of gluon degrees of freedom in the deconfined phase, 0(N 2 ). The fact that this already 
occurs for temperatures below (albeit close to) T c is likely due to residual finite-volume artifacts 
of the lattice simulations, which become particularly severe for second-order (or very weak first- 
order) phase transitions such as those of SU(2), SU(3) and SU(4). 

From the data for the trace of the energy- momentum tensor, it is then straightforward to 
obtain the other bulk thermodynamic quantities p, e and s by numerical integration: the results 
are shown in fig. [21 where the left, central, and right panels respectively show the pressure, the 
energy density and the entropy density (in units of the appropriate power of the temperature). 

Let us now discuss the continuum extrapolation of our lattice results. Since our numerical data 
are obtained from simulations of the pure-glue sector with the Wilson action eq. ([8]) , leading-order 
discretization effects are proportional to a~ 2 , i.e. to N^ 2 . To estimate the quantitative impact 
of deviations with respect to the continuum limit, we repeated our simulations of the SU(2) and 
SU(4) gauge theories on lattices with N t = 6, 7 and 8, keeping the temperature and physical 
volume fixed. Going from Nt = 6 to Nt = 8 is expected to reduce the dominating cutoff effects 
by a factor around one half, while keeping the computational costs at a constant signal-to-noise 
ratio limited. In fig. [3] we compare the values for A/T 3 for these two groups, as obtained from 
the three different sets of simulations: in the confined phase, discretization effects are very small, 
and compatible with the statistical errorbars. For this reason, one can safely assume that the 
systematic discretization effects affecting our Nt = 6 results (as well as finite-volume effects) are 
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Trace of the energy-momentum tensor and glueball gas 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



1.5 



0.5 



SU(2) 
SU(3) 
SU(4) 
SU(5) 
SU(6) 

all glueballs below the two-particle threshold 
contribution from the lightest glueball 




3.75 



0.85 



T/T 



Figure 1: Trace of the energy-momentum tensor (in units of T 3 ) as a function of T/T c , for the 
SU(./V) gauge groups studied in this work. The results displayed are obtained from simulations 
on lattices with Nt = 6 sites in the Euclidean time direction, except for the SU(3) Yang-Mills 
theory and for the points at the nine lowest temperatures of the SU(2) theory, obtained from 
simulations on lattices with Nt = 8. The dashed curve is the theoretical prediction for A/T 3 , 
assuming that the system can be described as a gas of non-interacting glueballs, and the dotted 
curve represents the contribution from the lightest state in the spectrum. 
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0.08 1 2i 2 




0.9 0.95 1 0.75 0.8 0.85 0.9 0.95 1 0.75 0.8 0.85 

tit t / r r/r 



Figure 2: The pressure (left panel) and the energy density (central panel), in units of T 3 , and 
the entropy density (right panel) in units of T 2 , as a function of T/T c , for the theories studied in 
this work. 

negligible with respect to the statistical errors and systematic uncertainties related to the scale 
determination. 

4 Comparing with a glueball gas 

Since the models that we are studying are pure gauge theories, the only physical states in the 
confined phase are massive glueballs. Hence, as a first approximation, it is reasonable to expect 
that the behavior of the thermodynamic observables shown in fig.[T]and in fig.[2]could be described 
in terms of a free relativistic gas of these glueball states. It is far from obvious that these glueballs 
should behave as free particles (at least for small values of N), and testing this assumption is the 
first goal of our analysis. As we shall see, this also requires an Ansatz for the glueball spectrum at 
high energies (in the vicinity of a possible Hagedorn-like transition), which, in turn, will allow us 
to discuss some subtle features of the phenomenological models used to describe this spectrum. 
Testing these string-inspired models is the second goal of our analysis. 

We performed the comparison of our data to the ideal glueball gas predictions in three steps: 

1. firstly, we assumed the gas to be dominated by the lightest glueball only; 

2. then, we included all the glueballs below the two-particle threshold, using the very precise 
numerical estimates available in the lattice literature; 

3. finally, we compared our data with the whole glueball spectrum, assuming a spectral density 
Ansatz inspired by the effective bosonic string model. 

A similar approach has been followed in ref. |16| for the SU(3) Yang-Mills theory in D = 3 + 1 
dimensions. 

In D = 2 + 1 spacetime dimensions, the pressure associated with a free, relativistic particle 
species of mass m is 
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Figure 3: Cutoff dependence of our results for the trace of the energy-momentum tensor (in units 
of T 3 ). These panels show the results obtained at the same temperatures, from simulations at 
three different values of the spacing a, corresponding to lattices with Nt = 6 (circles), 7 (squares) 
and 8 (triangles), for the SU(2) (left panel) and SU(4) gauge groups: the discretization effects 
appear to be comparable with or smaller than the statistical precision of our data. 

from which the other equilibrium thermodynamic quantities can be derived. In particular, the 
trace of the energy-momentum tensor can be written as 

^-^£i-p(-*?)~^-«-') <*> 

k=l 

(see also eq. (1A.8|) in the appendix |A|) . 

For the first two steps in the comparison of our lattice data to the glueball spectrum, we used 
the numerical values of the glueball masses and the parametrizations of the appropriate scaling 
functions, which are reported in ref. |14| . 

The curves in fig. [1] show the expected behavior of A/T 3 , for a gas of non-interacting glueballs: 
in particular, the dashed line is obtained summing the contributions from all glueball species 
(below the elastic scattering threshold) which are known from lattice spectroscopy calculations, 
while the dotted curve represents the leading contribution given by the lightest glueball. As 
already mentioned above, it is easy to see that both these curves fail to reproduce the data for 
T/T c larger than (approximately) 0.9. This has also been observed in recent, high-precision lattice 
computations of the equation of state for SU(3) Yang- Mills theory in D = 3+1 dimensions [16pi9j . 

Another feature, which is immediately manifest from the data, is the large separation between 
the bands of data corresponding to the SU(2) and the SU(A?" > 3) gauge groups. Since the value 
of the lowest glueballs (in units of T c ) is almost the same for SU(2) and for the other SU(AT > 3) 
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gauge groups, this gap must be the consequence of the fact that in the two cases the theories have 
different spectra, due to the aforementioned absence of C = — 1 states in the SU(2) Yang-Mills 
theory. This is an important difference with respect to the spectrum of the theories based on 
the other SU(iV > 3) gauge groups, admitting both C = +1 and C = — 1 states (not mutually 
degenerate). The fact that our SU(2) results for A/T 3 start to strongly deviate from those of the 
other groups at T/T c ~ 0.9 indicates that in this region the thermodynamics is likely dominated 
by effects due to the density of glueball states, rather than by just the lightest state in the 
spectrum. 

To describe the full glueball spectrum, various phenomenological models have been proposed 
in the literature: these include, in particular, bag-type models [37] and string-inspired models [9]. 
In the following we focus on the latter, and summarize their main features; besides the original 
work, the interested readers can find a discussion of its more recent generalizations in ref. [8]. 

In the original proposal by Isgur and Paton [9], glueballs are modelled as "rings of glue", i.e. 
as closed tubes of chromoelectric flux, which are described as closed bosonic string states. In 
particular, this implies that each glueball state corresponds to a given phonon configuration (i.e. 
to a given bosonic closed string state), and for each phonon combination there exists an infinite 
tower of radially excited states of increasing mass. This model can then be generalized, by 
including possible /c-glueball states (for N > 4), which correspond to closed fc-strings, metastable 
"adjoint string glueballs" (which become stable in the large- iV limit and may explain the splitting 
between the sectors of opposite C), as well as the finite thickness of the flux tube, which is usually 
modelled introducing an additional phenomenological parameter. This generalized version of the 
Isgur-Paton model turns out to be in remarkably good agreement with the low-lying spectrum of 
Yang- Mills theories, as calculated from first principles by means of lattice simulations [8]. 

An interesting feature of this model is that, essentially, these extensions lead to copies of 
the original spectrum, which are shifted towards higher values of the masses: thus, the thermo- 
dynamic contribution of the corresponding states is exponentially suppressed, except in a close 
neighborhood of a Hagedorn-like temperature. Furthermore, the correction to the spectrum due 
to the finite thickness of the flux tube becomes negligible for heavy glueballs, so that, as a first 
approximation, the glueball spectral density can be modelled in terms of the spectrum of a closed 
bosonic string (see the appendix [B] for details): 



For SU(iV > 3), the model predicts a further twofold degeneracy, accounting for the two possible 
orientations of the flux tube. 

Using this expression in eq. (|15p . assuming TW = T c and a Nambu-Goto string, for which 
T c 2 = 3a /it in D = 2 + 1 spacetime dimensions p one obtains the prediction shown in fig. HI 
In particular, the solid curve, which accounts for states of opposite charge-conjugation quantum 

2 It is interesting to note that, in D = 2 + 1 dimensions, the effective Nambu-Goto string model for con- 
finement predicts a numerical value for the ratio of the deconfinement temperature over the square root of the 
zero-temperature string tension T c j^fa approximately equal to one, in good numerical agreement with recent, 
accurate lattice determinations [27 , 38 . 




(16) 
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Trace of the energy-momentum tensor and string model 
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Figure 4: Same as in fig. [H but comparing our simulation results to the theoretical prediction 
including the contribution obtained from a bosonic string model for the glueball spectral density, 
as discussed in the text. The solid curve is the prediction for the SU(iV > 2) theories, while the 
dash-dotted curve accounts for the lack of C = — 1 states in the SU(2) gauge theory. 
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number, is relevant for the trace anomaly in SU(iV > 2) theories, while the dash-dotted curve, 
obtained including only C = +1 states, is expected to provide a more adequate description for the 
SU(2) theory. As the figure shows, both curves are in remarkably good agreement with our data, 
for all temperatures up to the region where the results corresponding to the different gauge groups 
start splitting from each other. This agreement is a strong piece of numerical evidence supporting 
the Isgur-Paton model, and, more generally, bosonic string models as effective theories for the 
confining regime of non-Abelian gauge theories. These results also support the hypothesis that 
glueball interactions are weak, and that the confined phase thermodynamics can be accurately 
approximated in terms of a relativistic gas of free massive bosons. 

In order to have a better feeling of the quality of the agreement between lattice data and the 
effective string prediction, it is also useful to compare our simulation results with the prediction 
that one would obtain, using the density of states of the open (rather than closed) string. This 
is displayed in fig. [51 where we compare the two curves to our data for SU(4) (the gauge group 
for which we performed the finest temperature scan) in the region where the contribution from 
the density of states of heavy glueballs dominates over the lightest ones: the precision of our 
simulations is sufficient to show that a model based only on open strings is clearly incompatible 
with the lattice results. However, our work does not rule out the possibility of modelling the 
glueballs in terms of a combination of closed and open string states in the adjoint (or in a higher) 
representation, as suggested in ref. [8]. 

Finally, in the close vicinity of T c , our results show that the contribution from heavier glueballs 
(or from interactions) becomes more and more important: this drives the change of behavior 
observed in the figures. It is interesting to note that, while the original Isgur-Paton model 
predicts exactly the same glueball spectrum for any number of colors N (except for the missing 
C = —1 states for N = 2), the extension discussed in ref. [8] predicts a dependence on N, related 
to the larger number of /c-glueball states which become available when N is increased. 

5 Conclusions 

In this work, we presented high-precision lattice results for the equation of state of SU(iV) Yang- 
Mills theories in 2 + 1 dimensions. We focused onto the confining phase, where the thermody- 
namics of these strongly coupled theories is expected to be described in terms of color-singlet 
hadronic states (glueballs). 

At low enough temperatures, the equilibrium thermodynamic properties are described well 
by a gas of non- interacting glueballs, with masses compatible with the results obtained from the 
accurate lattice determinations available in the literature. 

Close to the deconfinement temperature, however, this very simple model fails to reproduce 
the lattice data, and the contribution due to heavier glueball states has to be taken into account. 
The latter can be evaluated using a simple bosonic string model, like the one originally proposed in 
ref. [9] for the D = 3+1 case, or a refinement thereof [8]. The resulting equation of state, assuming 
that the Hagedorn temperature can be identified with the deconfinement temperature, and taking 
the effective string to be described by the Nambu-Goto model, is in very good agreement with 
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Closed versus open string model 




0.935 0.94 0.945 0.95 0.955 0.96 0.965 0.97 0.975 0.98 0.985 0.99 0.995 1 

T/T 



Figure 5: Dependence of the theoretical prediction for the equation of state on the string spectrum 
details: our simulation results for the SU(4) Yang-Mills theory are compared with the bosonic 
string prediction for the equation of state, assuming the glueballs to be modelled either as closed 
(solid curve) or open (dotted curve) strings. 



15 



the data from our lattice simulations. This gives further support to the validity of bosonic string 
models as effective theories for the confining phase in non-Abelian gauge theories. 

Our findings can be compared with those obtained in a similar computation for the SU(3) 
Yang-Mills theory in D = 3 + 1 dimensions [16], which also reported excellent agreement with 
an equation of state obtained extending the sum over known glueball masses with an exponential 
Hagedorn-like spectrum. One difference with respect to the latter work, however, is that in the 
present work we did not fit the value of the Hagedorn temperature Th'- in the D = 2 + 1 setup 
the deconfinement phase transition is a second-order one also for the SU(3) gauge theory — and 
a very weakly first-order one for SU(4) — and for continuous phase transitions it is expected that 
Th should be equal to T c . 
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A Ideal relativistic Bose gas in d + 1 spacetime dimensions 

The logarithm of the canonical partition function Z(T, V) of an ideal relativistic Bose gas is: 




where m is the mass of the boson and = 27r d / 2 /T(d/2) is the d-dimensional solid angle. 
Integration by parts yields: 





(A.2) 



where we set coshu 




1 + 2-%, and used the following integral representation 




(A.3) 
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for the modified Bessel function of the second kind of index v. 
In the thermodynamic limit the pressure is 



T 



( m 



2 x oo 



V 2ir 



d+l 

T \ ~ 



km J 



m 
f 



(A.4) 



The other equilibrium thermodynamics observables can be obtained from the above expressions 
for the pressure. For instance, the entropy density s is given by 



dp 
dT 



Similarly, the internal energy density e reads: 

T 2 8 



e = VdT lnZ = ~ P + sT 



(A.5) 



(A.6) 



Combining eq. (IA.4j) with the expression for the trace of the energy-momentum tensor in d spatial 
dimensions, = e — d ■ p, and using the recurrence relations of modified Bessel functions, one 
finds that the ideal Bose gas enjoys a remarkable identity: 



m 
2^ 



2 x ^ - 



E 

fc=i 



T 

km 



(A.7) 



namely, the trace of the energy-momentum tensor for the Bose gas in d spatial dimensions is 
proportional to the pressure Pd-2 of a Bose gas in d — 2 spatial dimensions: 



m 

&d = 7T~ Pd-2 

Z7T 

Finally, note that using the asymptotic expansion 

4i/ 2 -l 



valid for large |z|, one obtains: 



1 + 
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(A. 



(A.9) 
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(A.10) 



which, for d = 2, reduces to eq. ((T 
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B Spectral density of closed bosonic strings 



We assume that the glueballs can be modelled as "rings of glue", which are described, in the 
limit of large masses, by the Nambu-Goto model of closed bosonic strings. The mass spectrum 
in D = 2 + 1 spacetime dimensions reads: 

m 2 = Attg ^n L + n R - with: n L = n R = n, (B.l) 

where a is the string tension, the integers ul and n R describe the total contribution of the left- 
and right-moving phonons along the closed string, and the —1/12 term arises from the zero-point 
energy contribution. The degeneracy of these single-particle states is given by the number of 
partitions of til and n R (see, for instance, ref. [39]), hence the total degeneracy p(n) for the 
physical states is 

p{n) = -rr(n L ) n(n R ) = 7r(n) 2 , (B.2) 

where ir(n) denotes the number of partitions of n, and can be calculated using the generating 
function 

oo 1 oo 

IIt^^E-w 1 - ( b - 3 ) 

k=l y n=0 

For n large, one can resort to the Ramanujan asymptotic formula: 

.H^^exp^yf j , (B.4) 

which gives: 

In D spacetime dimensions, the Hagedorn temperature Th [6] is related to the string tension by: 



thus, for large m, 



so in D = 2 + 1 one gets: 



The spectral density as a function of the mass p(m) is defined via 

p(m) dm = p(n) dn . (B.9) 
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Using eq. (IB~7l) . one gets: 

dn= ? 3mdm 2 , (B.10) 



thus in D = 2 + 1 dimensions one obtains: 



3 



* m > e "" T "- (B11) 

The generalization to arbitrary D = d + 1 is straightforward: for a closed string in D spacetime 
dimensions, one finds: 

/OD (n) = 12( J D-2) D (^) e m / T ". (B.12) 
Combining this expression with eq. (1B~1I and eq. flBJOi one obtains: 

Mm) ^J^(^) D -\^. (B , 3) 
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